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We studied possible modifications of tlie neutrino - nucleon scattering rate due to the nucleon - 
nucleon collisions in the hot dense matter which we find in the supernova core. We show that the 
finite width of the nucleon spectral function induced by the nucleon collisions leads to broadening of 
the dynamical spin structure function of the nucleon, resulting in the reduction of the rate of neutrino 
- nucleon scattering via the axial vector current and making the energy exchange between neutrinos 
and nucleons easier. The reduction rate is relatively large ~0.6 even at density ~ lO^'^g/cm'^ and 
could have a significant impact on the dynamics of the collapse-driven supernova as well as the 
. . . cooling of the proto neutron star. 
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I. INTRODUCTION 



The collapse-driven supernova is supposed to be an outcome of gravitational collapse of a massive star ^ 8M0 at 
the end of its evolution. Although more than 60 years have passed since the original idea by Baade and Zwicky , 
^ ■ and the detection of neutrinos from SN1987A HQ confirmed that our scenario is correct on the whole, we have not 
10 \ yet figured out how this phenomenon occurs. 

■ These days many researchers of the supernova consider that neutrinos diffusing out of the proto neutron star play 
' a crucial role in heating the material behind the shock wave and expelling the outer layer of the star [^j6| . It turned 
out, however, that this mechanism is quite sensitive to the neutrino luminosity In fact, Janka and Miiller ||] 

showed by numerical simulations that only about 20% of increase in neutrino luminosity could lead to a successful 
Qs^ , explosion of a model which otherwise failed to explode. Hence mechanisms to boost the neutrino luminosity have 
been quested. Although many authors have been devoted in the study of convection in the core and it has been 
shown that the convection does help neutrinos heat matter, it is still controversial if the convection alone is sufficient 
for successful explosion or not l§-|l2] . It is also found that the sophistication of the numerical treatment of neutrino 



O ^ transport could increase the neutrino luminosity although the quantitative assessment of the difference it makes in 



the realistic context remains to be done llSHl 



^ ' On the other hand, our knowledge of the neutrino reaction rates in the hot dense medium is rather poor. In fact, 
. . , even in the most elaborate simulations of supernovac it has been assumed that the reactions of neutrinos with nucleons 
^ ■ are the same as in vacuum and the effect of surrounding matter is ignored p^-pj[|. However, the wavelength of a 
k>( \ 30MeV neutrino, for example, is longer than the mean separation of nucleons for the density ^ lO^^g/cm^, and the 
5_j ■ time scale corresponding to the same energy is roughly of the same order as the mean free time of nucleon between 
d collisions. Hence we have to study the spatial and temporal correlations of the matter, and it could be possible 
that the many body effects change the opacity for neutrinos considerably and we get the desired enhancement of 
neutrino luminosity and/or energy. In fact, efforts to find a possible modification of rates of neutrino reactions with 
nucleons have been made by several authors ||2C|-|3l|] . In these studies neutrino - nucleon scatterings are one of the 
targets, since it is one of the major sources of opacity for neutrinos. Reddy et al. [ p9pc| ] pointed out that taking 
a correct effective mass of nucleon into account changes the scattering rate in the high density regime ^ po which 
is the saturation density. Horowitz et al. |Q, Burrows et al. | |27| , p8[ |, Reddy et al. p9| , ^ and Yamada et al. ]3l| ] 
discussed the correlation effects due to the particle - hole excitation using a random phase approximation (RPA). In 
these studies the nucleons were assumed to be quasi-particles with the vanishing width of the spectral functions. On 
the other hand, Raffelt and his collaborators p4|-p6|| claimed that the neutrino scattering rate could be also reduced 
by losing the temporal correlations of the spins of nucleon, thus broadening the width of the structure function due 
to collisions of the nucleon with other nucleons surrounding it (see also p3|). This is a counter part for scattering 
of the so-called Landau-Pomeranchuk-Migdal effect [p2|-|34| for the bremsstrahlung. It was also pointed out that this 
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broadening of the structure function could enhance the inelasticity of the scattering and affect the energy spectra 
of neutrinos , since the width of the structure function is of the order of the matter temperature at a relatively 
low density ^ 10^'^g/cm'^ according to their estimate. In these studies they assumed that the structure function of 
nuclcon takes the Lorentzian form and inferred its width from the typical collision rate of nucleon with the overall 
normalization determined by the sum rule and the detailed balance imposed by hand. In addition, they only discussed 
long wave length limits. In this paper we study this coUisional effect on the neutrino - nucleon scattering rates on 
the field theoretical basis. The width of the spectral function of nucleon is evaluated by solving the Bethe-Salpeter 
equation in medium and calculating the imaginary part of the nucleon self-energy [p6H3q| , and then it is used to 
calculate the structure functions of nucleon as a function of the transfered four momentum In so doing the 

separable Yamaguchi potential |Q was used for simplicity. 

The paper is organized as follows. In the next section we formulate the basic equations. The results are presented 
in section 3. Discussion is given in the last section. 



II. FORMULATION 
A. scattering rates 

The scattering rates of neutrinos with nucleons are quite generally formulated as : 

R (El", E"""*, cos 6) = AGlEl''E°J'' [i?i(fc) (1 + cos 6*) + i?2(fc) (3 - cos 6i) 

- 2(i?™+i?°"*)i?5(fc)(l-cos0)] , (2.1) 

where and E°^* are the incident and outgoing neutrino energy, respectively, 9 is the scattering angle, and k is the 



four momentum transferred from neutrino to nucleon. In the right hand side of Eq. (2.1), the third term is usually 
much smaller than the other two, so that we ignore it in the following. The so-called dynamical structure functions 
Ri and i?2 are given in the non-relativistic limit which we assume in the following, as 



Ri{k)^ h'y j d\ e''^ {pN{x)pNm , (2.2) 
R2{k) ^ ^ Jd'x { s^^ix) siv(O) ) . (2.3) 

Here hy and Ha are the weak coupling constants of nucleons, ( ■ • ■ ) stands for the thermal ensemble average of the 
argument, and pn{x) and s^]^(x) are the density and the spin density of nucleon. Thus it is obvious that Ri (k) and 
i?2(fc) are nothing but their correlations in the matter. As is clear from the factor of the i?2(fc) in Eq- (P^), R2{k) 
is more important than Ri{k). Thus we will discuss in this paper the modification of R2{k) , or the spin-density 
correlation function, due to the nucleon - nucleon collisions in the supernova core. 



B. temporal spin-density correlation of a nucleon in random walk 

Just by the same argument as made by Knoll et al. ]39| ] for the momentum correlation, the spin-density correlation 
of the nucleon which is successively scattered by other nucleons with a mean collision rate, F, changing the direction 
of its spin via spin-dependent interactions, is given by ignoring the spatial non-uniformity as 

°° iFrl" 

S^\t) = {s\r)s\Q)) = e-l^^l5]^(.j„.^+„)„ . (2.4) 

Tl = 

In the above equation, n is the number of the scatterings during the time interval of t. Sm is the spin-density at the 
m-th collision, and (• • •)„ means taking the average of the argument with respect to the ensemble of the sequence of 
n scatterings. Taking the Fourier transform of S'^^{t), we obtain 

_ \^ / ^ k 2F"i?e(F + iuj)"+^ 

n— ^ ^ 

It was pointed out in the paper by Knoll et al. [ p9[ that the above summation of the Lorentz functions correspond to 
the summation of the Fcynman diagrams depicted in Fig. ^ in the closed time path formalism. 
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FIG. 1. The Feynman diagrams contributing to the spin density correlation function. 

If we further assume that the spin-density auto-correlation decreases by a constant rate, a, after each scattering, 
(Sm Sm+n)m = st)m, then we obtain 

y^'M- , .^^'^J ^'m^'j^n (2.6) 

r' = (l-a)r . (2.7) 

This implies that the higher order terms account for the difference between the collision rate, T, and the relaxation 
time of spin-density, T'. If the spin f lip i s fast enough, a « 0, these two rates become almost identical. In this case, 
the spin correlation function of Eq. (|2.8| ) is approximated by the first diagram in Fig. ^ Note that the propagator 
in that diagram is not a free propagator but a dressed one in the medium. The self-energy included in the dressed 
propagator comes from the nucleon collisions with other nucleons. In this paper this self-energy is evaluated by solving 
the Bethe-Salpeter equation for two nucleons as drawn schematically in Fig. || [|36|-^. In so doing, we assumed a 
separable Yamaguchi potential |^] in order to facilitate solving the Bethe-Salpeter equation. Although this is an 
oversimplification of the nuclear interactions, this is not a bad approximation as long as the low density regimes, 
Pb ^ 10^^ g/cm^^ is concerned, which is the case in this study, and gives some insight into what the collisional effects 
are like and how important they could be. Since the scatterings are suppressed due to the Fermi blocking in the high 
density regime, it is expected the scattering effects are most important in the low density regime. In addition, the 
resulting broadening of the structure functions of nucleon will affect the neutrino energy spectra which are formed in 
this low density regime . 
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FIG. 2. The Bethe-Salpeter equation of the T— matrix for two nucleons (upper) and the Dyson equation for the Green 
function of nucleon (lower). 



Some remarks are necessary. The above approximation is not applicable to the vector current contribution to the 
neutrino scattering rates, since it does not take into account the conservation law, and, as a result, the long wave 
length limit is not properly reproduced. In the case of the axial vector current, the spin density is not conserved if 



3 



one takes into account the nuclear tensor force. What is more important, the nuclear axial vector current is nearly 
isovector. Thus unless the matter is entirely composed of neutrons or protons, the axial vector component of the 
nuclear weak current is not a conserved one as long as the nuclear force is spin-dependent. 



C. approximations 



Here we summarize the formulation outlined above. The spin correlation functions are given by the dressed Green 
function as 



i?2(fc) 



hi 2 



G>ip + k) G<{p) 



(2.8) 



where the nucleon Green functions are defined as 



iG>{p) = P{p)[l ~ fN{p)] (2.9) 
iG<{p) - ^p{p)fN{p) ■ (2.10) 

fN{p) is the Fermi-Dirac distribution of nucleon. p{p) is the spectral function of nucleon given by the self-energy S(p) 



as 



ImT,{p) 



2Mn 



ReY.{p) 



1 ^ 

+ 


ImY,{p) 
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(2.11) 



Mn is the nucleon mass. The quasi-particle approximation, which is usually assumed, is obtained as a limit /ml](p) — > 
0. In the Hartree or Hartree-Fock approximation, the self-energy is real and the nucleon can be treated as a quasi- 
particle. Beyond those approximations, the self-energy has an imaginary part in general. It is obvious that the 
imaginary part of the self-energy is a width of the spectral function of the transfer energy Hence the spectral 
function becomes a S function in the quasi-particle approximation. 

In this paper, following Aim et al. we evaluate the imaginary part of the nucleon self-energy from the two- 
particle T-matrix which is given as a solution of the Bethe-Salpeter equation in the ladder approximation. 



T{pi,P2,Pi',P2',z) = 



V{pi,P2,Pl',P2') 

d^P3 d^p4 d^pz 



d^Pe 



(27r)3 (27r)3 (27r)3 (27r)= 



V{pi,P2,P3,P4) G2{p3,P4,P5,Pe,z) T{p5,pe,pi, ,p2' , z) 



(2.12) 



Here is defined as the product of two single-particle Green functions and calculated in the quasi-particle approxi- 
mation as 



G2{pi,P2,p'i,p'2,z) 



1 - /jv(Pi) - /iv(P2) 
z - E{pi) ~ E{p2) 



SHP1-P'1)S'{P2-P'2) 



(2.13) 



E{p) is the on-shell energy determined from the real part of the denominator of Eq. ( [2.11 ) and approximated as 
explained later. Eq. ( 2.12] ) is solved for the separable Yamaguchi potential, Va{p,p') — Wa{p)\aWa{p'), with Wa{p) — 

■ ■ - ■ ■ -12.3178 (MeV bn-^f/'^ for a = % and 



where the coupling constant and the effective range are 
Act = —14.6988 (MeV fm~^)^/^ for a = ^Si and 7 = 1.4488 fm^^, respectively. Then the imaginary part of the 



self-energy is given by 
/mS(p, uj + ie) 



d^p' 
(2^ 



[fN{E{p')) + gB{E{p')+uj)] ImTUP,p',P,p',E{p')+uj 



le) 



(2.14) 



Here again the quasi-particle approximation is used for a single-particle Green function. gsiE) = {exp[(£' — 2fi)/T] — 
is a Bosc-Einstein distribution function with the chemical potential of the nucleon pair. The subscript ex denotes 
the inclusion of the exchange term. The real part of the self-energy is also derived from the T-matrix in a similar way. 
Since it depends on the on-shell energy which is in turn determined by the real part of the self-energy, iterations are 
required to obtain the self-consistent value of the on-shell energy. For simplification we abandon this consistency in this 
paper and set the on-shell energy as E( p) — /2Mlj -f J7, where the effective mass and the potential U are taken 
from the relativistic mean field theory WlU42] . Using this on-shell energy and the corresponding imaginary part of the 
self-energy calculated by Eq. ( |2.14| ), we approximate Eq. ( ^.llj ) as p{p) = Imll{p)/{[p° - E{p)Y + [ImY.{p)/2]^}. 
Finally we evaluate Eq. (p^) with this approximate spectral function of nucleon. 
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III. RESULTS 



In Fig. ^, we show the imaginary part of the neutron self-energy as a function of the transfer momentum and 
energy for the density pb = 3 x lO^^g/cm'^ and temperature T — lOMeV and electron fraction — 0.4. The on-shell 
momenta roughly correspond to the ridge of the hill in this figure. The imaginary part of the self-energy is essentially 
a scattering rate of a neutron with another nucleon. It is evident from the figure that this is energy- and momentum- 
dependent and becomes as large as the average kinetic energy of the nucleon ^ 3T for small momentum transfers. 
/mS decreases as one goes to higher energies along the on-shell, since the scattering rate becomes smaller for higher 
energy incident nucleons. The self-energy of proton is almost the same as that of neutron in this case. For more 
asymmetric matter with smaller Y^, ImT,p becomes larger while /toI]„ gets smaller as seen in Fig. ^. This is simply 
because the scattering rate for the neutron - proton pair is larger than that for the neutron - neutron or the proton 
- proton pair. The greater degeneracy for neutron also contributes to the reduction of the neutron self-energy. The 
density-, temperature- and Yg-dependences of ImT, can be understood essentially in the same way. Since we solved 
the Bethe-Salpeter equation, the contribution from the deuteron pole, if any, is automatically taken into account. In 
fact, the parameter values of the Yamaguchi potential used in this paper are chosen so that the deuteron binding 
energy is reproduced in the vacuum. For the above cases, we confirmed no bound state can exist due to the medium 
effect as shown by Aim et al. For the lower density pb = 10^'^g/cm'^, however, the deuteron pole does exist and 
give tiny contributions to the imaginary part of the self-energy as shown in Fig. ^, where only the contribution from 
the deuteron pole was shown. The existence of the bound state is, thus, not relevant for the density, temperature and 
Ye of our interest. 

The spectral function is shown in Fig. ^ as a function of the transfer energy for different transfer momenta The on- 
shell momentum is located again at the peak of the spectrum. As mentioned above, this function would be a (S-function 
of the transfer energy if there were no interactions between nucleons. With collisions of two nucleons the spectral 
function is broadened, as expected, and the width of the spectral function is roughly the scattering rate for the given 
transfer momentum. It is seen that the width of the spectral function or the scattering rate of nucleon is dependent 
on the transfer momentum and the spectral function becomes more narrow-peaked as the transfer momentum gets 
greater, reflecting the decrease of scattering rates. It should be again emphasized that the width of the spectral 
function is of the same order for the low momentum transfer as the average kinetic energy of nucleon. Hence it is 
expected that the kinetics of neutrino - nucleon scattering will be also affected. The density-, temperature- and 
yg-dependences of the spectral function can be inferred from that of the imaginary part of the nucleon self-energy. 
As the density increases, the scattering rate of nucleon with another one is increased, thus leading to larger widths of 
the spectral function. Note that the greater degeneracy of nucleon counteracts via the Fermi blocking to reduce the 
scattering rate. The scattering rate is determined by this competition. As the temperature increases, the blocking 
effect is diminished. This enhances the scattering rate. However, the larger kinetic energy of nucleon lowers the 
scattering rates. Again the width is determined by the balance of these factors. As the matter becomes neutron rich, 
the width of the spectral function for neutron is decreased partly because the scattering rate for neutron - neutron 
pair is smaller than that for neutron - proton and partly because the neutron becomes more degenerate and the Fermi 
blocking reduces the scattering rate more effectively. 

Using the spectral functions obtained above, we calculated the nucleon structure functions with Eq. (^.^). The 
typical result is shown in Fig. |^ for R2{k, lu) as a function of the transfer energy uj for three different transfer momenta 
k. In this figure, three cases are compared. The thin dotted lines represent the structure functions for the free 
neutron, while the thin solid li nes p resent the case for which only the effective mass of nucleon was taken into account 
in the spectral functio n, En . (2.11). The smaller effective mass lowers the peak of the spectral function as pointed 
out by some authors although this effect is small for this low density {M^ = 890AIeV). It also tends to 

broaden the structure functions, since the recoil of the nucleon becomes larger as the effective mass of nucleon gets 
smaller. Again this is negligible in the low density regime considered here. On the other hand, the structure function 
is considerably broadened when the collisions among nucleons were taken into account, as shown with the thick solid 
lines in Fig. |^. The width is essentially given by the scattering rate of nucleons and taken over from the corresponding 
spectral function. The widths are not different for this low momentum transfers, which is qualitatively different from 
the results for the free nucleons and those for the quasi-particle approximation. As the density increases, the width 
becomes larger in general. As the temperature decreases, the part of the negative energy transfer is reduced, since the 
extraction of energy from the medium becomes more and more difficult as the nucleons become more degenerate. This 
is nothing but the detailed balance relation R{k) = e'^" R{—k) which is automatically satisfied in the formulation 
used in this paper. The width of the structure function becomes smaller for neutron and greater for proton as Y,, is 
decreased, as expected from the above discussions for the spectral function. 
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FIG. 3. The imaginary part of the self-energy for neutron as a function of the transfer energy E and the absolute value of 
the transfer three momentum. The baryonic density is pt = 3 x lO^^g/cm^, the temperature is T = lOMeV and the electron 
fraction is Ye = 0.4. 



The total scattering rate of the neutrino with the incident energy El" is given by the integration of the structure 
function with respect to the transferred momentum k and energy lo: 

R*°\E"') = / : i?(£;r, E°^\ COS 6) f 1 - /^(K™*)] 

(27r)3 io i-fc 2El"2E°-^* \ u , v ^ ! i jv\ v )\-. \ > 

with uf"°-^ = rnin{k, 2E™ — fc). The scattering angle 9 and the energy of the outgoing neutrino £'°"* are functions 
of E^J^ and k. The collisions of nucleons reduce in two ways the above scattering rates from the quasi-particle 
approximation. Not only the amplitude of the structure function is lowered, but the broadening of the structure 
function tends to put some fraction of the str ucture function outside the integration region. Ignoring the Fermi 
blocking of the outgoing neutrino in Eq. (^T|), we integrated the axial vector part for E^J^ = 3T and 

compared the result with the collision taken into account and the result for the free nucleon. For the model shown 
in Fig. ^ for example, we found that the scattering rate is suppressed by the factor of 0.41 for neutron and 0.38 
for proton, respectively. Note that the contribution of the effective mass is only a few percent in this case. For 
Pb = lO^^g/cm^, T = lOMeV and Ye = 0.4, we still found the reduction of - 0.6. 

The larger width of the structure function implies that the neutrino can exchange energy with nucleons more effi- 
ciently. In the supernova simulations performed so far, the neutrino ~ nucleon scattering was commonly approximated 
as the iso-energetic scattering, which means that this process does not contribute to the thermalization of the neutrino 
spectra. The above results stand by the claim by Raffelt and his company p^- |2^ , p5[ | that this may not be the case. In 
fact, the mean square root of energy transfer increased from ^ 5MeV to ~ 15MeV for the model of Fig. |^. Although 
this is still not a large value, it should be born in my mind that this process has about ten times larger cross section 
than the scattering with electrons. Hence it is expected the energy spectra of neutrinos, particularly muon- and 
tau-neutrinos, will be affected by this effect. Thus we need to somehow implement this effect in a numerical code and 
to study quantitatively the difference it will make in the realistic settings. 
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FIG. 4. The imaginary part of the self-energy for neutron (left) and proton (right) as a function of the transfer energy 
E and the absolute value of the transfer three momentum. The baryonic density is pb = 3 x lO^^g/cm^, the temperature is 
T = lOMeV and the electron fraction is Ye = 0.1. 




FIG. 5. The dcutcron contribution to the imaginary part of the sclf-cnorgy for neutron as a function of the transfer energy 
E and the absolute value of the transfer three momentum. The baryonic density is pb = lO^^g/cm^, the temperature is 
T = lOMeV and the electron fraction is Ye = 0.4. 
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IV. CONCLUSION 



We studied in this paper the possible effects of coUisions among nucleons on the reaction rates of the neutrino - 
nucleon scatterings. By solving the Bethe-Salpeter equation for two nucleons, we estimated the imaginary parts of the 
nucleon self-energy, which in turn give the spectral function of the nucleon in medium. Using this spectral function, 
we calculated the structure functions of nucleon, which are directly related to the neutrino - nucleon scattering rates. 
In the limit of no collision, the spectral function reduces to the delta function of the transfer energy and the structure 
function gives the reaction rates which have been commonly used in the literature thus far. It was shown that the 
collisional broadening of the structure functions might be the main factor in the low density regime considered in this 
paper to modify not only the neutrino reaction rates with the nucleon but also the neutrino spectra. 

It is obvious that the further study is needed on this issue. Although we used Yamaguchi potential to facilitate 
calculations, this is clearly an oversimplification of the dynamics of two nucleons. The self-consistency was also 
sacrificed for simplicity. Thus our results should be regarded as of qualitative nature. More i mpor tantly, we have to 
assess the corrections of the terms we ignored in estimating the the structure functions, Eq. (2.5). The in-medium 
vertex corrections should be also taken into account and are inferred from the Fermi liquid theory, for example . 
Though we neglected in this paper these corrections in order to make clear the effect of the width of the spectral 
function of nucleon, they could be as important and should be studied on the same basis in the future work. It should 
be also mentioned that the formulation in this paper is not conserving | |3l] , |45| , ^6[ | and cannot be applied to the vector 
current, whose long wave length limit is dictated by the baryon number conservation. The accomodation of this limit 
should be another subject of the future work. 

At present the neutrino transport is thought to play a crucial role in the collapse-driven supernovae. It is, thus, 
important to study the microphysics involved and, at the same time, to investigate its consequences to the supernova 
simulations, which we are planning to do with the numerical code recently developed by us to solve the Boltzmann 
equations for neutrinos fill. 




FIG. 6. The spectral function for neutron as a function of the transfer energy E for various absolute values of the transfer 
three momentum. The baryonic density is p;, = 3 x 10^"^g/cm^, the temperature is T = lOMeV and the electron fraction is 
Y, = 0.4. 
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FIG. 7. The dynamical structure function i?2(fc, w) for neutron as a function of the transfer energy u) for various absolute 
values of the transfer three momentum, k. The baryonic density is pi, = 3 x lO^'^g/cm'', the temperature is T = lOMeV and the 
electron fraction is = 0.4. The thin dotted lines represent the structure functions for the free neutron, while the thi n soli d 
lines present the case for which only the effective mass of nucleon was taken into account in the spectral function, Eq. (2.11). 
The thick solid lines show the results with the collisional effects taken into account for k = 20, 30, 50MeV/c from left to right. 
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